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ABSTRACT 


A  digital  computer  program  for  calculating 
either  continuous  ascent  or  stop-type  decompres¬ 
sion  schedules  is  described,  and  examples  of  appli¬ 
cations  are  given  The  formulas  used  for  continuous 
ascent  were  obtained  analytically  as  solutions  of 
differential  equations  relating  the  inert  gas  tension 
in  the  current  critical  tissue  to  the  safe  depth  and 
with  the  actual  depth  kept  equal  to  the  safe  depth 
at  all  times  after  an  initial  fast  rise  from  the  bottom 
to  the  safe  depth  Thus  the  rate  of  decompression 
of  the  critical  tissue  controls  the  rate  of  ascent  Gas 
tensions  in  nine  tissues  having  the  same  range  of 
gas  exchange  half-times  as  have  been  used  in  EDU 
reports  are  calculated  on  a  continuous  basis,  with 
the  one  having  the  deepest  safe  depth  being  the 
current  critical  tissue  The  stop-type  ascent  portion 
of  the  program  may  be  used  to  generate  a  staged 
ascent  using  the  same  parameters  for  comparison 
with  the  continuous  ascent  schedule  The  starting 
conditions  for  ascent  may  either  be  computed  by 
the  program  from  the  dive  history  or  be  communi¬ 
cated  to  it  as  a  subroutine  in  connection  with  an¬ 
other  program  The  program  may  be  used  either 
to  prescribe  ascent  schedules  or  to  analyze  dives 
for  which  the  history  is  known 
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COMPUTATION  OF  CONTINUOUS  DECOMPRESSION  SCHEDULES 

FOR  DEEP  SEA  DIVES 


INTRODUCTION 

Too  rapid  ascent  from  prolonged  deep  dives  can 
result  in  decompression  sickness  or  “the  bends,”  an 
illness  associated  with  inert  gas  bubble  formation 
in  the  blood  and  tissues  of  the  body  This  illness  is 
preventable  through  the  use  of  appropriate  decom¬ 
pression  schedules,  which  depend  principally  on 
the  dive  depth  and  its  duration  However,  many 
other  factors  that  may  affect  the  diver’s  physiology, 
such  as  exercise,  water  temperature,  and  the  effects 
of  breathing  different  gas  mixtures,  must  also  be 
considered  in  constructing  decompression  sched¬ 
ules  A  panoramic  view  of  recent  work  in  diving 
physiology  is  presented  in  a  book  edited  by 
Lambertsen  1 

Workman2  has  critically  reviewed  the  develop¬ 
ment  of  the  U  S  Navy’s  current  decompression 
schedules  3  These  are  based  on  the  principles  of 
inert  gas  exchange  as  developed  by  Haldane  and 
co-workers,4  and  on  their  determination  of  the 
amounts  by  which  the  inert  gas  tensions  in  the  tis¬ 
sues  may  safely  exceed  the  ambient  pressure  with¬ 
out  producing  the  bends 

The  currently  used  decompression  schedules 
involve  a  series  of  steps,  or  stages,  at  intervals  of  10 
ft  of  depth  The  diver  waits  at  each  stage  until  it  is 
safe  for  him  to  be  at  the  next  stage,  at  which  time 
he  ascends  to  it,  and  so  on,  until  the  surface  is 
reached  These  schedules  have  been  constructed 
to  minimize  the  decompression  time,  consistent 
with  preventing  the  bends,  and  have  been  proved 
safe  in  thousands  of  real  dives  However,  it  is  ap¬ 
parent  that  any  staged  decompression  schedule  is 
inherently  less  than  optimal,  because  as  soon  as 
the  diver  has  been  at  a  given  stage  for  even  a  short 
time,  it  would  be  safe  for  him  to  be  at  some  shal¬ 
lower  depth,  which  in  turn  would  give  an  increase 
in  the  outgassing  pressure  gradient  and  a  shorten¬ 
ing  of  the  time  required  for  further  decompression 
At  the  deeper  stages,  particularly  in  dives  with 
short  bottom  times,  some  tissues  may  actually  still 
be  absorbing  rather  than  releasing  inert  gases,  thus 


obligating  longer  stops  at  subsequent  stages  The 
importance  of  this  effect  increases  with  the  dive 
depth,  but  decreases  with  increasing  bottom  time 
The  potential  advantages  of  continuous,  rather 
than  staged,  ascent  schedules  have  previously  been 
recognized,  for  example,  by  Workman2  and  Brad- 
ner  and  Mackay  5  Shreiner  and  Kelly6  have  de¬ 
veloped  continuous  decompression  formulas  as¬ 
suming  constant  excess  saturation  pressure,  both 
for  constant  partial  pressures  of  oxygen  and  for 
constant  percentage  of  oxygen  in  the  breathing 
mixture 

The  purposes  of  this  report  are  (1)  to  present  the 
derivation  of  a  set  of  formulas  obtained  as  solutions 
of  differential  equations  that  include  the  assump¬ 
tion  that  the  depth  is  continuously  changing  in 
such  a  way  as  to  keep  the  diver  exactly  at  his  cur¬ 
rent  safe  depth,  and  (2)  to  describe  a  digital  com¬ 
puter  program,  simdive,  that  has  been  developed 
to  use  the  above  method  in  the  development  of 
more  nearly  optimal  decompression  schedules 
In  this  derivation,  we  do  not  consider  the  elab¬ 
orate  mechanisms  involved  in  bubble  formation, 
but  work  from  the  same  assumptions  used  in  de¬ 
veloping  the  standard  tables,  to  the  effect  that,  for 
any  particular  inert  gas  partial  pressure  in  a  tissue, 
there  is  some  corresponding  minimal  external  pres¬ 
sure  that  is  sufficient  to  prevent  bubble  formation  2 
The  depth  associated  with  this  minimal  pressure  is 
designated  the  safe  depth 

During  the  decompression  process  the  gases  are 
eliminated  from  the  body  and  the  required  ex¬ 
ternal  pressure  decreases  The  formulas  developed 
for  the  simdive  program  combine  the  equations 
for  the  rates  of  uptake  and  release  of  inert  gases 
with  equations  relating  the  internal  partial  pres¬ 
sures  to  the  safe  depth  A  schedule  is  produced 
such  that,  after  an  initial  rapid  ascent  to  the  safe 
depth,  the  actual  depth  is  kept  equal  to  the  instan¬ 
taneous  value  of  the  safe  depth 

The  formulas  developed  by  the  above  methods 
are  used  in  a  computer  subroutine  called  ascent 
Ascent  is  designed  to  be  called  by  simdive  or  an- 
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other  main  program  such  as  Moeller’s  standive7 
that  supplies  the  parameters  and  constants  needed 
to  express  the  conditions  of  the  dive  and  the  diver’s 
status  with  respect  to  his  tissue  inert  gas  partial 
pressures  at  the  beginning  of  the  ascent  Another 
subroutine,  descent,  provides  an  alternative 
method  for  generating  the  ascent  starting  condi¬ 
tions  from  the  dive  history 

The  simdive  program  and  associated  subrou¬ 
tines  were  originally  written  in  fortran  II  for 
compatibility  among  the  Honeywell-800  at  the 
SMC  and  the  IBM-7094  and  the  CDC-6600  com¬ 
puters  at  BNL  However,  the  version  presented  in 
this  report  includes  some  fortran  IV  features 
adapted  principally  for  use  with  the  CDC-6600 
computer 

MATHEMATICAL  BASIS 

It  is  assumed  that  the  driving  force  for  gas  ex¬ 
change  is  the  partial  pressure  gradient  between  the 
breathing  gas  mixture  and  the  gas  tension  in  the 
tissues  For  each  tissue  a  rate-limiting  gas  exchange 
factor,  A,  may  be  assigned  This  factor  lumps  the 
effects  of  several  parameters  It  depends  primarily 
upon  the  blood  perfusion  rate  per  gram  of  tissue 
and  the  solubility  of  the  gas  under  consideration 
in  the  tissue  relative  to  that  in  the  blood,  but  also 
reflects  variations  in  diffusion  rates  for  the  gases 
and  the  effects  of  temperature  changes  The  A’s  are 
related  to  tissue  gas  exchange  half-times  through 
the  relationship  XT  =  log,. 2  The  various  tissues  are 
regarded  as  being  in  parallel  for  gas  exchange  con¬ 
siderations,  and  thus  are  treated  independently 
For  consistency  with  other  work2’7  nine  tissues 
having  a  spectrum  of  A’s  corresponding  to  gas  ex¬ 
change  half-times  ranging  from  5  to  240  min  are 
used  These  half-times  were  arbitrarily  chosen  to 
give  the  desired  spread  of  values,  and  do  not  neces¬ 
sarily  correspond  to  any  particular  anatomical 
tissues  or  organs  The  factors  relating  tissue  gas 
tensions  to  the  safe  depth  are  also  functions  of  the 
half-times  Because  of  the  different  rates  of  uptake 
and  release  of  the  gases  among  the  tissues,  there 
will  usually  be  one  for  which  the  safe  depth  at  a 
given  moment  is  greatest,  and  this  will  be  called 
the  critical  tissue  If  two  or  more  tissues  tie  for  this 
distinction,  the  one  with  the  longer  half-time  is 
taken 

The  plan  for  ascent  is  first  to  ascend  as  rapidly 
as  is  practical  (60  ft/nun)  until  the  actual  depth 
equals  the  safe  depth  for  the  critical  tissue  From 


this  time  on,  the  actual  depth  is  kept  equal  to  the 
safe  depth  for  the  current  critical  tissue  Control 
passes  from  one  tissue  to  the  next  when  the  curves 
representing  their  individual  safe  depths  cross 
One  of  the  programming  problems,  as  discussed 
below,  is  to  predict  the  time  for  changing  the  criti¬ 
cal  tissue  For  dives  involving  short  bottom  times, 
the  rapidly  exchanging  tissues  are  the  controlling 
ones  at  the  beginning  of  ascent,  and  control  passes 
to  progressively  slower  tissues  as  the  surface  is  ap¬ 
proached  For  very  prolonged  dives,  the  slower  tis¬ 
sues  may  be  rate-controlling  from  the  beginning  of 
ascent 

For  comparison,  a  variation  of  the  program  pro¬ 
duces  a  staged  ascent  In  this  variation  the  diver 
remains  at  a  given  stop  until  a  time  such  that  by 
ascending  to  the  next  even  10-ft  level,  at  the  maxi¬ 
mum  rate,  he  will  arrive  at  the  new  depth  at  the 
instant  this  becomes  his  safe  depth  This  is  slightly 
different  from  the  standard  tables,  in  which  the 
ascent  time  is  not  considered,  but  gives  a  fairer 
comparison  with  the  continuous  method 

The  notation  to  be  used  is  given  in  Table  1 

Insofar  as  possible,  the  derivations  that  follow 
are  m  terms  of  a  single  tissue,  and,  for  simplicity  of 
expression,  subscripts  are  omitted  when  possible 
without  ambiguity  In  the  computer  program  the 
subscripts  are  needed  to  distinguish  among  the 
various  tissues  and  gases 


General 


If  (P—  Q)  is  the  driving  force  and  A  is  the  rate- 
limiting  factor  for  inert  gas  exchange,  then 

(1) 

or,  in  standard  linear  equation  Form, 

§+Ae=w, 

(2) 

also 

AD=AD0-DD=ADo-  f  ‘  Iidl , 

A  o 

(3) 

P=G(AD  +  33)  =  G(D-  J  lRdt ) 

(4) 

Safe  Depth 

Workman2  mentions  several  formulas  that  have 
been  used  to  calculate  safe  inert  gas  partial  pres¬ 
sures  in  tissues  for  given  depths  In  the  present  no- 
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Table  1 

Symbol 

Explanation 

Units 

a  (  jj ) 

7  th  gas  tension  in  J th  tissue  (range 

for/,  1-3,  for  J,  1-9) 

ft  of  water 

Ci(U) 

I th  gas  tension  in  critical  tissue,  L 
(7=1  is  used  for  total  inert  gas 
tension,  1  =  2  for  nitrogen,  and 

f  =  3  for  helium) 

ft  of  water 

Q.UJ)  o 

<2  (./,/)  at  start  of  current  step 

ft  of  water 

a  uJh 

Safe  value  for  (7  (JJ )  at  current 

depth 

ft  of  water 

Pin 

Partial  pressure  of  7  th  gas  in 

breathing  mixture 

ft  of  water 

G(I) 

Fraction  of  7th  gas  in  breathing 

mixture 

— 

t 

Time  from  start  of  current  step 

min 

R 

Rate  of  change  of  depth,  taken  as 

positive  for  ascent 

ft/min 

AD 

Actual  depth  (gauge) 

ft  of  water 

SD(J ) 

Safe  depth  for  /  th  tissue  (gauge) 

ft  of  water 

A  D0 

Depth  at  start  of  current  step 

(gauge) 

ft  of  water 

DD 

ADo  —  AD=  f  Rdt,  change  of 

depth  for  step 

ft  of  water 

D 

AD( i  +  33,  total  external  pressure 

ft  of  water 

A  U) 

Turnover  rate  constant  for  J  th 

tissue,  A  =  (loge  2)/half-time 

inin_1 

AK{J) 

Q.  (JJ )«  at  surface 

ft  of  water 

AC{J) 

A  A  K(  J  )/ft  of  depth 

— 

Note  All  depths  and  pressures  are  measured  in  feel  of 
sea  water,  absolute,  unless  gauge  depth  is  speci¬ 
fied  (For  fresh  water  dives  D  —  ADq  +  34,  and  the 
AK’ s  and  dC’s  are  multiplied  by  1  025,  the  den¬ 
sity  of  sea  water  ) 


or, 

SD=SD0--Q2Z£1  (8) 

Ascent  Formulas 

As  mentioned  above,  the  ascent  schedule  has 
two  phases,  first  a  rapid  ascent  to  the  minimum 
safe  depth,  then  a  slower  ascent  maintaining  the 
current  safe  depth  The  former  involves  a  constant 
rate  of  change  of  depth,  whereas  the  latter  involves 
a  continuously  changing  rate  The  formulas  for 
constant  rate  are  simpler  and  are  derived  first 
These  foimulas  are  also  used  for  descent  and  for 
computing  staged  ascent  schedules 

Constant  Rate 

With  R  constant, 

f  Rdt  —  Rt  ,  (9) 

-'o 

Therefore,  substitution  of  Rt  into  Eq  (4)  gives 

P=  G(D  —  Rt) ,  (10) 

and  Eq  (2)  becomes 

^+A  d=\G(D-Rt),  (11) 

for  which  the  solution  is* 

Q.  =  exp(  —  A/)|  J  AG(D  —  Rt)  exp(Xt)  dt  +  C  j 

(12) 


tation,  the  following  formula  is  equivalent  to  the 
one  used  for  constructing  the  tables  of  recom¬ 
mended  M  values  in  ref  (2) 

Qt  =  AK+AC-AD  (5) 

By  rearrangement  of  Eq  (5),  the  safe  depth  for 
a  given  value  of  Q.  is 

SD=(Q_ -AK)/AC  (6) 

Thus,  if  Q  changes,  the  *SD  changes  proportionally 


SDq  —  SD 


(Qo-AK) 

AC 


( Q.-AK ) 
AC 


Integration  and  simplification,  with  C=Q  0,  yield 
Q  =  GD[  1  —  exp(  — A03—  1  +e-xp(-A<)]  + 

A 

Q.0exp(  — A<)  (13) 

In  this  form  it  is  seen  that  the  third  term  on  the 
right  represents  the  function  Q.  for  G=0,  l  e  ,  if  the 
diver  breathes  pure  oxygen  Similarly,  the  first 
teim  represents  the  build-up  in  tissue  inert  gas 
partial  pressure  that  would  occur  if  77  =  0  and 
Q  o  =  0  The  middle  teim  corrects  for  the  effect  of 


(Qo-Q) 

AC 


(7) 


*The  symbolism  exp  ( X )  is  used  here  to  represent  ext  where  e 
is  the  base  of  the  natural  logarithms,  2  718 
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the  rate  of  ascent  For  computing  purposes,  a  re¬ 
arrangement  of  Eq  (13),  with  the  exponential 
terms  collected,  is  used 

Q_=GD-  ^±z}l  +(  Q  o  -  CD-  exp(  -  At) 

(14) 

For  constant  depth,  with  R  =  0,  Eq  ( 1 4)  becomes 

Q=,CD+(Q0-CD)exp(-\t)  (15) 

The  time  required  for  Q_  to  reach  some  other 
value,  Q_  t,  at  constant  depth  may  be  calculated  by 
solving  Eq  (15)  fort 

^loge[(£o  -CD)/(Q_t-GD)} 

A 

provided  that  GD  ^£7i<Q.oor£7o<(7*<  GD 


The  first  t0  is  obtained  by  substituting  the 
approximation 

exp  (  — At)=r  1  +At—  \  (At)2  (22) 

into  Eq  (19),  and  the  result  into  Eq  (18),  and 
solving  the  resulting  quadratic  expression  for  t  In 
subroutine  ascent,  the  critical  tissue  at  the  time  of 
intersection  is  determined  by  selecting  the  tissue 
for  which  the  substitution  involving  Eq  (22)  gives 
the  least  value  for  t,  then  the  exact  value  of  t  is  cal¬ 
culated  for  the  critical  tissue  only,  by  using  Eq 
(21)  When  the  diver’s  actual  depth  equals  his  safe 
depth,  another  set  of  equations  involving  the  use  of 
a  variable  ascent  rate  is  needed  The  derivations 
for  these  equations  follow 

Variable  Ascent  Rate 

To  maintain  AD  =  SD(L), 


Intersection  of  Fast  Ascent  and  Safe  Depth  Curves 

To  find  the  time  of  intersection  of  the  fast  ascent 
and  safe  depth  curves,  we  set 


(23) 


and 


AD=SD  (17) 

and  substitute  Eq  (3)  and  Eq  (8) 

AD0-Rt  =  SD0-22_lJll  (i8) 


But,  from  Eq  (13), 

(Qo— 2)  =  (2o-C£)[l-exp(  — A0]  + 

GR  [Ai  —  1  +  exp(  —  At)]  ^ 

A 

Substitution  of  the  right-hand  member  of  Eq 
(19)  into  Eq  (18)  eliminates  Q.  and  gives  an  ex¬ 
pression  involving  only  t  and  variables  whose 
values  are  known  at  the  start  of  a  given  step  Al¬ 
though  this  expression  cannot  be  solved  analyti¬ 
cally  for  t,  a  value  of  t  may  be  obtained  to  any 
desired  accuracy  through  successive  iterations  of 
Newton’s  aproximation  formula  8 
If  it  is  desired  to  solve  F(t)  =  0  for  t. 


DD 


=  flRdt  = 


[£(L,l)o-£(L,l)] 
AC(L) 

Substitution  of  Eq  (24)  into  Eq  (4)  gives 


/*=G  D  — 
and  Eq  (2)  becomes 

mjj) 


=XjG[ 


dt 


D- 


A  C(L) 


+*jQU,n 


i- 


[£(/.,l)o-<?QL,l)J 

AC(L ) 


(24) 


(25) 


(26) 


In  order  to  obtain  an  integrable  form  of  Eq 
(26),  it  is  necessary  first  to  solve  for  Q(L,  1)  With 
the  understanding  that  only  the  total  inert  gas 
tension  in  the  critical  tissue,  £7  (L,l),  is  under  con¬ 
sideration,  the  notation  may  be  simplified 


t — Iq 


m 
no  ’ 


(20) 


where  to  is  an  initial  guess  for  /,  and  F'(t)  =  dF(t)/dl 
The  approximation  obtained  by  the  first  solution 
of  Eq  (20)  for  t  is  used  as  the  to  for  the  next  itera¬ 
tion,  etc  ,  until  the  answer  converges  within  the 
desired  limits  Thus  t  is  obtained  from  successive 
iterations  of 


f-+Aa=AC(°-!r)+^f’  (27) 

or,  in  standard  form 

J  §-+(‘-£)x<i=xc(d-j?)  <28> 


AD0 -SD0  -  Rt  +  {( Q_  o  -  GD)  [  1  -  exp(  —  A;)]  +  GR{Xt  -  1  +  exp(  -  A/)]/A}A4  C 
—  R  +  {X(Q_o  —  GD)cxp(—\t)  +  GR[[  —  exp(  —  Xt)]}/AC 


(21) 
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In  Eq  (28),  let 

K  /,  g\  AC~G 
Al_l  AC/  ~  AC  ’ 

k*=g(d-jZ)’ 

„  *2  G(D-AC-Qo) 

Ks~  Kx~  AC-C 


Substitution  of  the  right-hand  member  of  Eq  (36) 
into  Eq  (26)  gives 

=\jG,\d-  tQ-  (2:,!)°  —  A~3]  [  1  —  exp(-XLA~iO)| 
t  AC(L)  ) 

(37) 

From  the  definition  of  A3  it  is  readily  shown  that 


Then 


^L+\k1q=xk2 

(29) 

Solving  Eq  (29)  gives 

Q,  =  exp(  —  XKxt)  I  f  XK2  exp(AATit)  dt  4-  C 
1-  Jo 

(30) 

Integrating  between  limits  gives 

r  1  r  , 

Q,  =  exp(  -  XKxt)  |  [exp(A Ail)  —  1  ]  +  C  j 

(31) 

For  t  =0,  C=  Q.  0, 

Q,  =  A’3[1-  exp( — A A"i<)]  +  Q.  0  exp( —XKil) 

(32) 

=  +  ( Q  0  —  A3)  exp(  —  A  Kxt) 

(33) 

(38) 

AC  Gx  V  ' 

With  this  substitution,  Eq  (37)  reduces  to 

dQ-(J,I)+^jdUJ) 

=A,Ci|  ^  +  {D~  -§)«P(-^10]  (39) 

with  the  solution 

a(7,/)=exp(-AJ0{  f\jGj[§-  + 

Gl  ° G-\  ^  exp(- - j  exp(X  J  ^dt  +  Cj  (40) 

Integration  of  Eq  (40)  gives 


for  (l  (L,  1),  1  e  ,  the  total  inert  gas  tension  in  the 
critical  tissue  [Eq  (33)  is  not  directly  applicable 
for  other  tissues  or  for  individual  inert  gases  ] 

From  Eqs  (3),  (24),  and  (33),  we  obtain 

Z)Z)  =  (Q,0  — /f3)[l  —  exp(— \A'1t)]//4C  (34) 

To  solve  Eq  (34)  for  t,  let 

K4  =  (Q0-JCa)/AC , 


Q.UJ)—  exp(-AJo{-^p(^[exp(AjO-1]  + 

GXD-K ,  v  ) 

\j  —  A£/A‘1(exP[(^-/— —  U  )  +  (41) 

au,n= -  exp(-Aj<)]+ 

X  fexp(  —  XLKxt)  —  exp(  — Aj-t)]  j  + 


then 

t=  -  loge(l -DD/KJ/XKt  (35) 

Thus  at  this  point,  for  the  critical  tissue,  L,  we 
have  analytical  solutions  for  Q_  ,  R,  and  DD  as 
functions  of  time,  and  for  time  as  a  function  of  the 
change  in  depth,  DD,  all  with  AD  —SD  Whether 
Eq  (34)  or  Eq  (35)  is  to  be  used  depends  on 
whether  t  or  DD  is  known,  respectively  We  pro¬ 
ceed  to  solve  Eq  (26)  for  the  other  tissues  and  for 
the  component  inert  gases 
From  Eq  (33), 

£(L,l)o  — £(Z.,1) 

=  [£(^,1)o-A3][1-  cxp(-Ajr10]  (36) 


Cexp(  —  Xjt)  (42) 

for  t  =  0,  C=Q_U,I) 0 
Finally 


2(y,/)=g-[A 3+ 


A  j{GxD  —  K3) 
Xj—XlKi 


exp(  —  A^Ait) 


+ 


(43) 

An  indeterminancy  arises  in  Eq  (42)ifA  j—XiKx, 
with 


exp(  —  \iKxt)—  exp(  — A  jt)  _  0 

Xj—X^Ki  0 


(44) 
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no 


This  indeterminancy  may  be  removed  through  an 
application  of  L’HopitaPs  rule 

Lim  I  exp(  — A^AV)—  exp(-Ajt)~[ 

.\,Ai  I  \j  —  Ai,Ai  I 

—  l  exp(  —  Xl,K\()  (45) 

In  ascent,  a  modification  of  Eq  (43)  involving 
substitution  of  Eq  (45)  is  used  if 

lAy-A^IClO-e 

Otherwise  an  expression  equivalent  to  Eq  (43)  is 
used 


and 


t=zt0- 


m 

no 


(51) 


Iq  this  case  a  satisfactory  first  estimate  for  to  is  ob¬ 
tained  by  substituting  the  approximation 


exp(A")zr  1  +X 
into  Eq  (49)  and  solving  for  i 

[> SD(L)0-SD(J)o \ 


to  —  : 


Q.  (L)o  —  GiD  Q.(J)o-GiDx 

AL - - -  AJ 


(52) 


(53) 


AC(L) 


AC(J ) 


Change  of  Critical  Tissue 

Except  in  saturation  dives,  when  the  critical  tis¬ 
sue  is  the  one  with  the  longest  half-time,  the  criti¬ 
cal  tissue  changes  whenever  the  safe  depth  curve 
of  one  of  the  slower  tissues  crosses  that  of  the  cur¬ 
rent  critical  tissue  It  is  useful  to  be  able  to  predict 
when  this  will  occur,  and  this  may  be  done  through 
the  use  of  the  formulas  derived  below 

The  critical  tissue  changes  when,  for  J^>L  (as¬ 
suming  that  the  tissues  are  ordered  by  increasing 
half-times), 


SD(J)  =  SD(L) 


From  Eq  (8), 


SD(J)  o- 


=  SD(L)  o- 


SD(L)a-SD(J)o+ 


[Q.U)o-Q.U)} 

ACU ) 

re  Wo -£(/.)] 


(46) 


AC(L ) 

mj)o-nj)} 

AC{J) 


(47) 


rewo-ew] _ n 

AC(L) 


(48) 


Substitution  from  Eq  (36)  and  Eq  (43)  with 
Gi—  G\  gives 


SD(L)0-SDU)o  + 


The  same  expression  for  to  is  obtained  if  the  time  of 
intersection  of  the  tangents  at  t  =  0  for  the  two  curves, 
SD(J)  and  SD(L)  as  defined  by  Eq  8,  is  taken 

PROGRAM  DESCRIPTION 

Listings  and  logic  diagrams  for  program  simdive 
and  its  subroutines  are  given  in  Appendix  B 

As  the  “main”  program,  simdive  reads  in  the 
initializing  data  and  calls  the  appropriate  sub¬ 
routines  The  first  two  columns  of  the  first  card  of 
a  data  deck  determine  whether  the  run  will  begin 
with  descent  or  will  go  directly  to  ascent,  and 
whether  new  half-time,  AK,  AC,  and  gas  mixture 
parameters  or  those  for  the  preceding  deck  are  to 
be  used  The  balance  of  this  card  through  column 
72  may  be  used  to  give  any  alphanumeric  infor¬ 
mation  desired  to  identify  the  dive  (For  consist¬ 
ency  with  other  reports,  the  half-times  for  nine 
tissues  usually  used  are  5,  10,  20,  40,  80,  120,  160, 
200,  and  240  mm  )  If  descent  is  used,  and  the 
diver  is  not  at  the  surface  when  the  end  card  is 
reached,  the  program  automatically  proceeds 
with  ascent  without  first  returning  to  simdive 
When  the  surface  has  been  attained  in  ascent, 
there  is  a  return  to  simdive,  which  proceeds  to  call 
table  and  plot,  then  starts  over  or  stops,  depend¬ 
ing  on  whether  another  dive  follows 


M(  J)o-*i]H  -  exp(-Xjt)]+  expt-X^Q]  [eWo_/r,1|,  _ 


mj) 


AC(L) 


=  0 


Again,  Eq  (49)  cannot  be  solved  explicitly  for  l,  but  a  numerical  solution  is  obtainable  by  the  same 
method  as  used  above  for  the  fast  intersection  If  F(t)  is  the  left  member  of  Eq  (49), 


mj) 


AC(L) 
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Subroutine  descent  is  used  primarily  to  com¬ 
pute  the  inert  gas  partial  pressures  in  the  tissues  at 
the  beginning  of  ascent  For  each  step  it  reads  in 
the  time  and  depth  for  the  end  of  the  step  and  the 
index  for  the'  gas  mixture  used  during  the  step  (In 
an  alternative  version,  the  actual  gas  mixture  frac¬ 
tions  are  read  in  here  )  The  starting  time  and 
depth  for  a  given  step  are  assumed  to  be  those  for 
the  end  of  the  immediately  prior  step  Subroutine 
tension  is  called  for  computing  the  inert  gas  par¬ 
tial  pressures  in  the  tissues  at  the  end  of  each  step 
Descent  may  also  be  used  to  compute  the  safe 
depths,  partial  pressures,  and  safe  partial  pressures 
for  the  decompression  steps  in  any  dive  in  which 
the  rates  of  change  of  depth  are  all  constants 
Finally,  descent  may  be  used  to  study  the  con¬ 
tinuation  of  the  decompression  process  after  the 
diver  has  surfaced,  for  example,  to  obtain  the  starting 
partial  pressures  in  the  tissues  for  repetitive  dives 

Ascent  may  be  used  as  an  independent  sub¬ 
routine,  or  may  be  joined  with  descent,  as  is  done 
in  the  listing  given  in  Appendix  B  Ascent  inserts  a 
dummy  first  step  of  zero  duration  to  print  out  the 
starting  conditions  for  ascent  It  automatically  de¬ 
termines  the  gas  mixture  to  use,  based  on  the  start¬ 
ing  depth  for  each  step  The  safe  depths  for  all 
nine  tissues  are  computed,  and  the  tissue  having 
the  greatest  safe  depth  is  selected  as  the  critical 
tissue  If  the  distance  from  the  actual  depth  to  the 
safe  depth  for  the  critical  tissue  is  >0  01  ft,  the 
program  calls  subroutine  fastriz,  which  takes  the 
diver  to  the  current  safe  depth  Otherwise  the  slow 
ascent  method  is  used,  in  which  the  actual  depth 
and  safe  depth  are  kept  equal  Normally  a  step  is 
an  interval  of  10  ft,  but  if  a  change  of  critical  tis¬ 
sues  or  a  gas  mixture  boundary  is  encountered 
first,  the  information  as  of  the  intervening  depth  is 
printed  out  If  the  rate  of  ascent  is  >  1  0  ft/min, 
ascent  proceeds  for  1  min  rather  than  10  ft  for  that 
step  In  the  printout  the  critical  tissue  for  the  step 
is  marked  with  an  asterisk  If  the  “stop”  is  due  to  a 
change  in  critical  tissue,  the  next  critical  tissue  is 
also  marked 

Subroutme  fastriz  finds  the  time  of  intersection 
of  the  actual  depth  curve  with  the  safe  depth 
curve,  assuming  ascent  at  the  maximum  rate, 
which  is  taken  as  60  ft/'min  unless  another  value 
is  read  in  by  simdive 

Subroutine  plot  produces  a  graph  of  the 
descf.nt-ascent  depth-time  curve  The  ordinate 
scale  is  based  on  maxdeep  A  caption  and  ordinate 
labels  for  0,  ’A,  l/z,  %,  and  1  times  maxdeep  are 


printed  out  The  time  scale  is  adjusted  to  keep  the 
plot  on  two  pages,  through  factors  entered  as  data 
statements  Once  the  scales  have  been  determined, 
plot  proceeds  by  treating  each  line  of  printout  as 
a  tune  interval,  and  examining  the  end  times  for  the 
dive  steps  in  sequence  If  a  step  end  time  falls  with¬ 
in  the  current  time  interval,  the  corresponding 
point  is  labeled  with  an  asterisk  When  a  step  is 
encountered  that  is  later  than  the  current  time  in¬ 
terval,  this  interval  is  printed  out,  the  labels  for  the 
steps  that  fell  within  the  interval  are  “erased”  (l  e  , 
labeled  blank),  and  the  program  proceeds  to  the 
next  time  interval  For  intervals  containing  points 
printed  out,  the  time  of  the  last  point  in  the  inter¬ 
val  is  printed  After  the  last  step  of  the  dive,  the 
ordinate  scales  and  captions  are  printed  out  again 
Subroutine  table  lists  the  starting  and  ending 
time  and  depth  and  the  gas  mixture  for  each  step 
The  times  are  converted  from  minutes  to  hours 
and  minutes  for  use  in  table 

Subroutine  tension  uses  the  constant  rate  for¬ 
mula,  Eq  (14),  for  computing  the  inert  gas  partial 
pressures  in  the  tissues 


Input 


INPUT-OUTPUT 


The  input  data  deck  is  described  in  Table  2 


Output 

Simdive  assigns  consecutive  serial  numbers  to 
the  dives  in  a  given  run  This  dive  number,  the 
identifical  ion  field  of  card  1 ,  and  the  date  of  the 
run  are  printed  at  the  beginning  of  the  ouput,  be¬ 
fore  table  and  before  plot  If  they  are  new  pa¬ 
rameters,  H,  AK,  and  AC  are  printed  For  descent, 
the  dive  number,  step  number,  prior  depth,  end- 
of-step  depth,  oxygen,  nitrogen,  and  helium  per¬ 
centages,  safe  tissue  pressures  and  ()  values  (total 
inert,  nitrogen,  and  helium)  for  each  of  nine  tis¬ 
sues,  the  step  time  interval,  the  cumulative  dive 
time,  the  rate  of  ascent,  and  the  safe  depth  for  the 
critical  tissue  are  printed  out  for  each  step  The 
critical  tissue  is  marked  with  an  asterisk  For 
ascent,  the  first  five  gas  mixtures  (oxygen,  total 
inert,  nitrogen,  and  helium)  and  the  depths  for 
which  each  is  to  be  used,  and  the  maximum  rate 
of  ascent  aie  printed  out  The  dive  number,  etc  , 
as  listed  for  descent  are  printed  except  that  the  (7 
values  are  followed  by  the  step  time  interval,  the 
cumulative  ascent  time,  the  cumulative  dive  time, 
and  the  average  ascent  rate  for  each  step,  instead 
of  the  items  listed  for  descent  Since  the  actual 
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Table  2 

Input  data  deck 

Entry 

Meaning 

Card  1  Col  1 

Blank  or  0 

Use  preceding  parameters 

1 

Read  new  parameters  (Cards  2-14) 

Col  2 

Blank,  0,  or  1 

Continuous  ascent  only 

2 

oescent,  then  continuous  ascent 

3 

descent,  then  staged  ascent 

Col  3-72 

Alphanumeric 

Dive  identification 

If  col  1,  card  1  blank  or  zero,  skip  to  card  15 

Card  2  8x  9F8  0 

Real  type  numbers 

H{J)J=  1,9 

Card  3  8x  9F8  0 

It  11  11 

AC(J  ),./=  1,9 

Card  4  8x  9F8  0 

11  (I  i? 

mj\j=  1,9 

Cards  5-14  Col  9-16 

M  If  II 

Nitrogen,  gmix(3,j) 

Col  17-24 

il  il  il 

Helium,  GMtx(4,j),y=  1,10 

Card  15  Col  17-24 

il  ll  il 

Initial  time,  t(1) 

Col  25-32 

n  ll  n 

Initial  depth,  depth(I) 

Col  33-37 

Integer 

maxoeep  (for  scaling  ordinate) 

Col  38-45 

Real  type  number 

rate  (maximum  ascent  rate) 

Cards  16-18  8x  9F8  4 

II  11  II 

Initial  tissue  pressures, 

0.  (IJ),  7=1,9,  7=1,3 

If  col  2,  card  1  blank,  0,  or 

1 ,  skip  to  card  M+  2 

Cards  19-7V*  Col  1-3 

Blank 

Col  4-8 

Real  type  number 

End  of  step  hour 

Col  9-13 

11  If  11 

End  of  step  minute 

Col  14-20 

II  II  11 

End  of  step  depth 

Col  46-50 

Integer  1-10 

Index  for  step  gmix 

Card  JV+ 1  Col  1-3 

END 

Indicates  end  of  descent  deck 

Card  N+2  Card  1  of  next  data  deck  or  end  of  file  card 

*N  may  be  any  number  19-217,  with  one  card  per  step  in  descent 

depth  and  the  safe  depth  are  kept  equal  after  the 
initial  fast  ascent,  the  safe  depth  is  not  printed 
separately 

After  the  diver  has  surfaced,  or  the  198th  step,  a 
tabulation  of  the  step  numbers,  step  times,  begin¬ 
ning  and  end-of-step  depths,  beginning  and  end- 
of-step  times,  and  the  gas  mixture  for  each  step  are 
printed  by  subroutine  table 

Subroutine  plot  prints  a  graph  of  depth  vs  time 
for  the  dive 

Execution  time  of  the  program  on  the  CDC- 
6600  averages  about  2  sec  per  dive  for  the  600-ft 
dives 

APPLICATIONS 

Figure  1  compares  two  ascent  schedules,  one  for 
staged  ascent  and  the  other  for  continuous  ascent, 
both  as  calculated  with  the  simdive  program  fol¬ 
lowing  descent  to  450  ft  for  1  hr  In  both  dives  it  is 
assumed  that  the  breathing  gas  mixture  is  air 


during  descent  to  75  ft,  then  96%  helium,  4% 
oxygen  for  the  rest  of  descent,  while  at  450  ft,  and 
during  ascent  to  200  ft,  after  which  air  is  again 
used  Although  there  is  very  little  difference  be¬ 
tween  the  two  schedules  for  the  early  stages,  the 
difference  gradually  accumulates,  and  the  diver 
surfaces  2  hr  sooner  by  the  continuous  method 
In  Figure  2  the  standard  table  schedule  for  a 
dive  to  190  ft  for  1  hr  on  air  is  compared  with  as¬ 
cent  schedules  calculated  by  the  simdive  staged 
and  continuous  methods  Although  the  simdive 
staged  schedule  moves  through  the  deeper  stages  a 
few  minutes  slower  than  does  the  standard  table, 
it  suggests  that  less  time  is  needed  at  the  10-ft  stage 
and  that  the  diver  could  surface  about  20  min 
sooner  Use  of  the  continuous  schedule  would  save 
another  20  min  in  this  dive 

In  Figure  3  a  schedule  for  ascent  as  computed 
by  the  simdive  method  is  compared  with  the 
schedule  given  in  the  diving  manual  table  for  ex- 
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Figure  1  Comparison  of  staged  and  continuous  ascent 
schedules,  both  as  calculated  by  simdive  for  a  helium  dive 
to  450  ft  for  60  min 


min  ,  2  3  4 

TIME  IN  HOURS 


Figure  3  Comparison  of  continuous  ascent  schedule  cal¬ 
culated  by  simdive  with  the  corresponding  schedule  from 
the  diving  manual  for  a  dive  to  300  ft  for  1 80  mm  on  air 


ns CIHT  GAS 
MIXTURES 


Figure  2  Comparison  of  staged  and  continuous  schedules 
calculated  by  simdive  with  the  staged  schedule  from  the 
standard  tables  for  a  dive  to  190  ft  for  60  mm  on  air 


Figure  4  Continuous  ascent  schedules  for 
mixed  gas  dives  to  600  ft 


treme  exposures3’11 103  for  decompression  following 
an  air  dive  to  300  ft  for  3  hr  Here  it  is  seen  that 
for  the  shallow  stages  the  safe  depth  as  calculated 
by  the  present  method,  and  assuming  that  the 
diving  manual  schedule  is  being  followed,  is  in 
close  agreement  with  the  depths  prescribed  in  the 
standard  tables  However,  there  is  a  large  differ¬ 
ence  at  the  deep  stages  Thus  the  difference  be¬ 
tween  the  two  schedules  is  due  more  to  the  differ¬ 
ent  relationship  assumed  between  tissue  tension 
and  the  safe  depth  for  the  deep  stages  than  to  the 
difference  between  staged  and  continuous  ascent 
schedules 

Figure  4  shows  a  family  of  continuous  ascent 
schedule  curves  for  dives  to  600  ft  for  various  bot¬ 


tom  times  ranging  from  6  min  to  12  hr  The  gas 
mixtures  used  are  tabulated  at  the  right-hand  side 
of  the  graph  Changes  in  the  gas  mixtures  account 
for  some  of  the  sharp  changes  in  slope,  while 
change  of  critical  tissue  accounts  for  others  As 
the  240-min  half-time  tissue  becomes  more  nearly 
saturated,  it  becomes  the  critical  tissue  earlier  in 
ascent  After  about  20  hr  bottom  time,  further  in¬ 
creases  in  bottom  time  do  not  increase  the  ascent 
times 

DISCUSSION 

The  above  derivations  provide  the  formulas 
needed  for  producing  any  desired  decompression 
schedule,  either  with  staged  ascent  or  with  a  con- 
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tinuous  ascent  in  which  the  actual  depth  is  kept 
equal  to  the  minimum  safe  depth  at  all  times  after 
an  initial  fast  ascent  to  the  safe  depth,  with  the 
assumption  that  the  safe  total  inert  gas  partial 
pressure  in  the  critical  tissue  is  linearly  related  to 
the  depth  according  to  Eq  (5)  If  some  other  re¬ 
lationship  is  required,  the  formulas  will  have  to  be 
modified  accordingly  Other  possible  improve¬ 
ments  include  increasing  the  number  of  tissues 
considered,  or  even  using  a  continuous  spectrum  of 
gas  exchange  half-times  Either  of  these  changes 
would  produce  a  somewhat  smoother  schedule 
than  does  the  present  method,  in  which  the  rate 
changes  abruptly  when  the  critical  tissue  changes 
The  schedules  produced  by  the  present  method 
are  optimal  in  the  sense  that  they  lead  to  the 
shortest  possible  decompression  schedules  consist¬ 
ent  with  the  assumptions  made  concerning  the  re¬ 
lationship  between  tissue  tension  and  safe  depth 
However,  until  any  differences  in  schedules  sug¬ 
gested  by  this  method  have  been  carefully  tested 
with  real  ascents,  they  canot  be  recommended  for 
use  in  practice 

The  practical  difficulties  involved  in  executing  a 
complicated  schedule  approximating  the  optimal 
one  is  recognized  Analogue  computer  methods9*10 
are  being  developed  that  may  eventually  automate 
the  controls  and  make  any  desired  schedule  feasible 
The  principal  value  of  the  continuous  ascent 
schedules  at  present  is  to  indicate  the  areas  in 
which  improvements  over  the  present  schedules 
may  or  may  not  be  expected  to  be  possible  The 
program  can  also  be  used  to  analyze  given  dive 
schedules  and  may  be  of  use  in  pinpointing  the 
reason  if  bends  result  during  an  actual  dive 

REFERENCES 

1  CJ  Lambertsen,  Editor,  Proceedings  of  the  Third  Sym¬ 
posium  on  Underwater  Physiology,  Williams  &  Wilkins 
Co  ,  Baltimore,  1967 

2  RD  Workman,  Calculation  of  Decompression  Schedules 
for  Nitrogen-Oxygen  and  Helium-Oxygen  Dives,  U  S  Navy 
Experimental  Diving  Unit  Research  Report  6-65, 1965 

3  US  Navy  Dicing  Manual,  Navships  250-538,  Washing¬ 
ton,  1963 

4  A  E  Boycott,  G  C  C  Damant,  and  J  S  Haldane, 
Prevention  of  compressed  air  illness,  J  Hygiene  8,  342- 
443  (1908) 

5  H  Bradner  and  R  S  Mackay,  Biophysical  limita¬ 
tions  on  deep  diving  Some  limiting  performance  ex¬ 
pectations,  Bull  Math  Biophys  25,251-72  (1963) 

6  H  R  Schreiner  and  P  L  Kelly,  Computation 
methods  for  decompression  from  deep  dives,  pp  275- 
99  in  Ref  1 


7  G  Moeller,  Stan  dive,  fortran  Solution  of  Decom¬ 
pression  Equations,  U  S  Naval  Submarine  Medical 
Center  Report  A65,  1966 

8  The  authors  are  indebted  to  Arthur  Harris,  Applied 
Math  Dept  ,  BNL,  for  suggesting  the  use  of  this 
method  for  obtaining  near-exact  approximations  in 
the  otherwise  intractable  equations  (19)  and  (49) 

9  BA  Hills,  A  thermal  analogue  for  the  optimal  de¬ 
compression  of  divers  Theory,  Phys  Med  Biol  12, 
437-44  (1967),  Construction  and  use,  Ibid  445-54 

10  RA  Stubbs  and  D  J  Kidd,  Computer  analogs  for 
decompression,  pp  300-1 1  in  Ref  1 


APPENDIX  A 

Principal  Variables  Used  in  SIMDIVE  and 
Subroutines,  in  Order  of  First  Use 


Input 

Variable  format  Explanation 


Program  SIMDIVE 

LDIVE 

Serial  number  of  dive  in  this  run 

K 

12 

Determines  type  of  run  (/T=12, 
descent  with  new  parameters, 
K  =  1 0,  11,  ascent  with  new  pa¬ 
rameters,  K=  2,  descent  with 
former  parameters,  K  =  0,  1,  as¬ 
cent  with  former  parameters  ) 

head(7) 

7A10 

Dive  identification 

h(9) 

9F8  0 

Half-times  in  minutes  for  nine 
tissues 

ac(9) 

9F8  0 

Slope  of  curve  relating  safe  tension 
to  depth 

ak(9) 

9F8  0 

Safe  total  inert  gas  partial  pressures 
in  tissues,  in  feet  of  sea  water 

dkgon(9) 

KJ)  =  ^e2/H(J) 

gmix(4,  10) 

2F8  2 

Fractional  partial  pressures  of  ni¬ 
trogen  and  helium  for  10  mixtures 
(total  inert  =  N2  +  He,  02  =  1  — 
total  inert) 

fcmix(3,  10) 

Ratio  of  individual  inert  gas  to 
total  inert  gas 

i(200) 

F8  2 

t(  1  )=  starting  time,  in  minutes 

depth(200) 

F8  2 

depth(I)  =starting  gauge  depth 
in  feet  of  sea  water 

MAXDEEP 

15 

Maximum  depth  to  be  shown  on 
graph 

RATE 

F8  0 

Maximum  rate  of  ascent 

q(9,  3) 

9F8  4 

Total  inert,  nitrogen,  and  helium 
partial  pressures  (absolute)  in  nine 
tissues,  in  feet  of  sea  water 

Input 

Variable  format 


Input 

Variable  format 


11 


Explanation 


Subroutine  descent 


VEND 

A3 

Signal  for  end  of  descent  data 

TllOUR 

F5  0 

Hours  component  of  step  time 

TMIN 

F4  2 

Minutes  component  of  step  time 

depth(200) 

F8  0 

Gauge  depth  for  step,  in  feet  of  sea 
water 

kc(200) 

15 

Index  number  identifying  gas  mix 
for  step 

dtime(200) 

Time  interval  for  step,  in  minutes 

saeeq(9) 

Safe  total  inert  gas  tension,  in  feet 
of  sea  water 

DD 

Change  in  depth  during  step,  in 
feet  of  sea  water 

rate2 

Average  rate  of  change  of  depth  for 
step,  in  ft/nun 

safed(9) 

Safe  depth  (gauge)  for  total  inert 
gas  tension  in  tissues 

Subroutine  ascent 


STIME 

Starting  time  of  ascent,  in  minutes 

CUTIME 

Accumulated  ascent  time 

BUM 

Next  depth  boundary  for  change  of 
gas  mixture 

DOD 

Absolute  external  pressure  (gauge 
depth  +  33),  in  feet  of  sea  water 

DB 

Distance  from  start  of  step  to  next 

BUM 

DEVEN 

Distance  from  start  of  step  to  next 
even  10-ft  level 

Explanation 

pti(9)  Provisional  time  of  intersection  of 

J th  tissue  with  critical  tissue  safe 
depths 

Subroutine  FASTRIZ 

All  variables  defined  above  or  by  program  statements 


Subroutine  PLOT 


INDEX 

Determines  ordinate  scale 

jy(5) 

Ordinate  labels 

ATIME 

Time  in  minutes  from  t(1)  to  end 
of  last  step 

TSCALE 

Scaling  factor  for  abscissa 

point(IOI) 

Ordinate  location  to  be  plotted 

ipoint(50) 

point’s  actually  used 

KLS 

Time  in  minutes  for  abscissal  loca¬ 
tion  K 

Subroutine  TABLE 

ntime(4) 

n  hme(  1 ),  (2),  beginning  step  time 
in  HR  MIN 

ntime(3),  (4),  end  step  time  in 

HR  MIN 

Subroutine  TENSION 

PPRES 

Inert  gas  partial  pressure  in  breath¬ 
ing  gas  mixture,  in  feet  of  sea 
water,  absolute 

PRDL 

Correction  for  effect  of  ascent  rate 
on  dQ/dt 

Listings  and  Logic  Diagrams  for  Program  simdive,  Subroutine  DESCENT,  Subroutine  ASCENT, 
Subroutine  fastriz,  Subroutine  plot,  Subroutine  table,  and  Subroutine  tension 
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A  digital  computer  program  for  calculating  either  continuous  ascent  or  stop  type 
decompression  schedules  is  described,  and  examples  of  applications  are  given.  The 
formulas  used  for  continuous  ascent  were  obtained  analytically  as  solutions  of  differ¬ 
ential  equations  relating  the  inert  gas  tension  in  the  current  critical  tissue  to  the 
safe  depth  and  with  the  actual  depth  kept  equal  to  the  safe  depth  at  all  times  after  an 
initial  fast  rise  from  the  bottom  to  the  safe  depth.  Thus  the  rate  of  decompression 
of  the  critical  tissue  controls  the  rate  of  ascent.  Gas  tensions  m  nine  tissues  hav¬ 
ing  the  same  range  of  gas  exchange  half-times  as  have  been  used  m  EDU  reports  are 
calculated  on  a  continuous  basis,  with  the  one  having  the  deepest  safe  depth  being  the 
current  critical  tissue.  The  stop  type  ascent  portion  of  the  program  may  be  used  to 
generate  a  staged  ascent  using  the  same  parameters  for  comparison  with  the  continuous 
ascent  schedule.  The  starting  conditions  for  ascent  may  either  be  computed  by  the 
program  from  the  dive  history  or  be  communicated  to  it  as  a  subroutine  in  connection 
with  another  program.  The  program  may  be  used  either  to  prescribe  ascent  schedules 
or  to  analyze  dives  for  which  the  history  is  known. 
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